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Following a relativistic heavy ion collision, the quark-gluon plasma produced eventually undergoes 
a chiral phase transition. We assume that during this phase transition one can describe the dynamics 
of the system by the linear a model and that the expansion can be thought of as mostly radial. 
Because the a model is an effective field theory there is an actual momentum cutoff (Landau pole) in 
the theory at around 1 GeV. Thus it is necessary to find ways of obtaining a covariantly conserved, 
, renormalized energy-momentum tensor when there is a cutoff present (which breaks covariance), 

• in order to identify the effective equation of state of this time evolving system. We show how to 

^— j I solve this technical problem and then determine the energy density and pressure of the system as 

a function of the proper time. We consider different initial conditions and search for instabilities 
which can lead to the formation of disoriented chiral condensates (DCCs). We find that the energy 
, density and pressure both decrease quickly, as is appropriate for a rapidly cooling system, and that 

■ the energy is numerically conserved. 

!>■ ■ 
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O . I. INTRODUCTION 

oo ■ 

' 

QQ ' During a relativistic heavy ion collision, it is possible to create a quark-gluon plasma, which will then cool and 
[ expand, leading to hadronization. At around the same time as hadronization, the system also undergoes a chiral 

■ phase transition which breaks the SU(2) x SU(2) symmetry and which then gives physical masses to the pion and 
0^ other particles. The out of equilibrium dynamics of this phase transition can be very interesting if the expansion is 

, fast enough so that the phase transition resembles a quench. In this case the effective pion mass can go negative for 
O i' short periods of proper time which leads to a temporary exponential growth of low momentum domains where the 
' isospin can point in a particular direction (DCC). This is a topic of current interest since DCCs may provide 
a signature of the chiral phase transition. In an earlier paper ^ the time evolution of this chiral phase transition, 
(-H I assuming a uniform spherically symmetric expansion into the vacuum, was considered. The 0(4) linear a model in 
• • I the large-A^ approximation, which incorporates both noncquilibrium and quantum effects, was studied. A range of 
initial conditions was analyzed to determine which initial states could lead to the formation of DCCs. This radial 
expansion is interesting because it is expected at late times in the plasma evolution, and also maximizes the cooling 
^ ' rate and therefore the possibility of the phase transition resembling a quench. 
5^ , In this paper we go beyond previous work |^ and study the proper time evolutions of the effective hydrodynamic 
collective variables of the system, namely the renormalized pressure and energy density. We want to see to what extent 
previous intuition coming from studying classical hydrodynamical models of particle production such as Landau's 
model is confirmed. This requires the solution of a new technical issue, namely how one obtains a covariantly 
conserved energy- momentum tensor when one has an actual cutoff in the theory. As discussed in earlier work |^,^ , 
the linear a model is an effective theory with a renormalized coupling constant of order 10, in order to agree with low 
energy pion properties. This leads to a maximum cutoff of around 1 GeV in order to avoid problems at scales of the 
order of the Landau pole. It has been shown earlier |3| that as long as the momentum cutoff determined from the 
sum over mode functions is below the Landau pole, there is a regime of cutoffs where the continuum renormalization 
group flow is obeyed so that one can safely say that one is in the continuum regime. For this class of problems, 
renormalization methods based on formal schemes such as dimensional regularization approaches are not very useful. 
Here we have a real cutoff in physical momentum as a result of using an effective field theory, and we have to be very 
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careful in order to obtain the correct covariantly conserved energy-momentum tensor which leads to our definitions 
of comoving energy density and pressure variables. The tools we use to correctly renormalize the energy-momentum 
tensor are first to introduce the physical cutoff A and then to analyze the divergences in all the components of the 
energy-momentum tensor using an adiabatic expansion of the mode functions. By comparing our results with a 
covariant point-splitting approach [Q, we can identify those terms which should survive in the covariant limit. This 
technique allows us to obtain finite energy densities and pressures which then enable us to study the effective equation 
of state for the evolving plasma. 

In this paper, we derive the energy-momentum tensor for the 0(4) linear a model, and show how to regularize and 
renormalize it in order to obtain the physical energy density and pressure. We study initial conditions which have 
instabilities that lead to DCC formation and also stable initial configurations, and examine the energy density and 
pressure in both situations. We also consider the conservation of the energy-momentum tensor as a check on our 
numerical methods. 

The paper is organized as follows. In Sec. l| we discuss the model and coordinate system used, and show how 
to construct the energy-momentum tensor. Then in Sec. Ill we explicitly describe the scheme used to renormalize 



this tensor in order to perform a numerical simulation with a physical momentum cutoff. In Sec. IV we numerically 
examine the conservation of energy and the equation of state of the system. Finally in Sec. ^ we discuss the results 
and provide some concluding remarks. 

II. ENERGY-MOMENTUM TENSOR 
A. Equations of Motion and Coordinate System 

We consider the 0(4) linear a model, with the classical action given by 



where the mesons are in a 0(4) vector representation 

$,-((7,7?). (2.2) 

The Lagrangian density is 

C['f^,x] = lg,A^WMx)WMx)] - lx{x)'S>U^) + yX(2^) + j^x'ix) • (2.3) 

In order to quantize the system, we construct the generating functional of the connected Green's functions and carry 
out the path integral to obtain the following quantum effective action (to leading order in N) ||] 

iN 

r[0„x;j.] = S[(P,,x;M + —TT\nGo\x,x';x) , (2.4) 

where 



Go\x, x'; x) = [□ + X{x)] S\x, x')/^/^^ , (2.5) 



and 



Mx) = (Mx)), 

X{x) = ixix)) . (2.6) 

By varying the effective action with respect to and x, we obtain the following equations of motion for the mean 
field (pi 

[a + xix)]M^)^Mx) = HS,o , (2.7) 
and the constraint equation for the effective mass squared x 

X{x) = A[-«2 + <^2(^) „ iNGoix, x; x)] . (2.8) 
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Notice that Gq{x,x';x) depends on x through its definition, and it fulfills the following dynamical equation 



[□ + xix)] Go{x, x'- x) = 5\x, x')/^^^). (2.9) 

The current has only a non- vanishing component in the zero (cr) direction in order to give mass to the pion . There 
are three parameters in the model: the mass of the pion ra = m,r; the value oi H = f-Kfn^, which gives the vacuum 
expectation value (tr) = f-n (using PCAC); and the coupling constant A, which is determined by fitting to low energy 
■n — IT scattering data |^,|| . 

The picture one gets from hydrodynamical simulations of heavy ion collisions is that the energy density is initially 
in a Lorentz contracted disk which expands first in the longitudinal direction and becomes three dimensional at late 
times. At early times one expects the velocity to scale approximately as z/t, where z is the longitudinal direction, since 
the effective longitudinal size goes to zero with center of mass energy. This leads to the energy density approximately 
becoming a function only of the longitudinal fluid proper time variable = (t^ — z^)^/^. At later times when the 
expansion is more spherical and the initial distribution looks more like a "point" , one expects that the velocity scales 
as V = r/t and the energy density then becomes a function of the spherical fluid variables: 

1 ^t + r\ 

77 = arctanh(r/t) = ^^^i j , (2-10) 

where t = r cosh 77 and r = rsinh?7 [^. We restrict the range of these variables to the future light cone, namely 
< T < cx) and < r/ < 00. The coordinates (r, 77, 0, </?) are useful to describe a spherically symmetric expansion of a 
plasma when one is in a hydrodynamical scaling regime. A spherical expansion provides the fastest possible expansion 
rate (and therefore cooling rate) of the quark-gluon plasma and thus enhances any nonequilibrium effects that are 
based on the idea of a rapid quench. Since complete inhomogeneous evolutions are at the edge of or beyond what is 
presently numerically possible, this spherically symmetric expansion provides the other extreme when compared with 
the slower, purely longitudinal expansion studied earlier 

In terms of this coordinate system, Minkowski's line element 

ds^ ^dt^-dr^- {dO^ + sin^ Od^^^ ) (2.11) 

is given by 

ds^ = dr^ - T'^{djf + sinh^ ?7 dO"^ + sinh^ r/sin^ dip^) , (2.12) 

from which we can read off the metric tensor 

g^i, = diag (1, -r^, -t^ sinh^ 77, -r^ sinh^ rj sin^ 9) . (2-13) 

The metric in this coordinate system is of the Robertson- Walker form: 

gi_,^ = diag (1, -a^(r), -a^(T) sinh^ ry, -a^(r) sinh^ ry sin^ 6*) , (2-14) 

which corresponds to a Ricci flat cosmological model with uniform expansion a(T) = r, and hyperbolic spatial sections, 
i.e., with curvature k = —1. 

It will be convenient to introduce the conformal time u given by 



a-\T') dr' , (2.15) 



or equivalently 



T 



e" 



m 



(2.16) 



where we choose m = ttItt, the only mass scale in the system. The metric in the coordinate system with conformal 
time (u, T], 9, ip) has the form 

g^^ = C{u) diag (1, -1, - sinh^ 77, - sinh^ 77 sin^ 9) , (2.17) 

with 
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C7(u) = a2(r) = r2=e2"/m2, 

y/^ ^ C'^{u)smh^r]sme . (2.18) 

In a hydrodynamical model, all the expectation values depend only on the proper time of the system. We therefore 
will assume that the mean fields (jji and x S'^e only functions of u, that is (f>i = 4>i{u) and x = xi"^)- We then write 
the full quantum field $i in terms of its expansion about its mean value, as follows: 

77, 6*, (j)) = + 77, 6*, 0), (2.19) 

where are the quantum fluctuations, which include both vacuum and thermal excitations. The equations of motion 
are then 

[□+X(w)]0.(a:) =0, (2.20) 
where x is the four vector x — (w, ry, 0, (p). Then for Go{x, a;'; we find 

Go{x,x';x) = * (Tc{0i(a;), 0j(a;')}) , 

where Tc corresponds to a M-ordered product, following the closed-time-path formalism of Schwinger ||6|. 

In order to solve the wave equation for the quantum fluctuations, we follow Parker and Fulling ff, and write a 
mode expansion for as follows 



(2.21) 



•^0 Ira 

with 
and 

A(3)X,„(f) + (s2 + l)y,irn{x) = , 

where A*^^^ is the Laplacian of the three dimensional hyperbolic spatial sections of curvature n — —1. Then the mode 
functions gs{u) satisfy the following differential equation 

g^iu) + [s^ + Ciu)xiu)] gsiu) = , (2.22) 

and for Go{x, x; x) we find Q 

{<fi)^-iGo{x,x;x)^C-\u) __(2n, + l)|5,(«)|2, (2.23) 

once we have chosen a particular vector state with respect to which we shall be taking expectation values. We choose 
an initial state such that the pair densities are zero, and the particle number density is finitely integrable with respect 
to the corresponding integration measure, namely 





slni) 




(5(s'- 


- s) Sii> Smm' > 




shri) 


= ("s + 




d{s' - s) Sii> S„im 




'm' ) 


= Ps % 


5{s' ~ 


s) Sw Smm' = , 




'ra' ) 


= P*s 


5{s' - 


S) Sw 5mm' = . 



Here we choose the initial particle number density to be a thermal distribution 
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where Wg(ito) = + C{uo)x{uo). Notice that we can choose the pair density to vanish {ps = 0), since one has the 
freedom to make a Bogohubov transformation at uq so that this always remains true. The boundary conditions on 
the mode functions gs{u) is that they correspond to the positive frequency adiabatic mode functions initially. We 
therefore choose 



9siuo) 

The following rescalings will be useful: 



gs{uo) = ^/^y2uJs{uQ), 



2ws(uo) 



Uu) = C~^'^{u)p,{u), 
x{u)^C-\u)x{u), 

In terms of the scaled variables, the equations of motion can be written 



with the gap equation being 



v\u) + pi{u) + N 



gs{u) = 0, 
Piiu) = j^{u) 



27r- 



-{2ns + l)\gs{ut 



(2.24) 

(2.25a) 
(2.25b) 
(2.25c) 
(2.25d) 

(2.26) 
(2.27) 

(2.28) 



B. Construction of the Stress-Energy Tensor 

The energy-momentum tensor is defined by 

T^, = (1 - 2e)(V^<I',) (V^D,) + [2^ - 5p..9"''(V„$,) (V^<i>,) 
- 2e$.(V^V,$,) + 1^5^,$, + ^(1 - 30x.9p.*' 

Equivalently, we can write this equation as 1^ 

T,,, - (V^$,) {V,<^,)+^{g^,0 - V^V,)$2 _ (£[$„x] + h<^i) , (2.30) 



where is given by Eq. ( |2.3| ). We take the expectation value of the energy-momentum tensor in the thermal 

initial state chosen, and make use of the mode expansion given in Eq. (2.21) to examine the components of T^^. 
After some algebra, we find for the T^u component 

C(ii)(r™)-^[p' + (x + i-i26p?] 

+ (64 - ^)PiPi-]tPi ^ ^ 



'(2n, + 1)<^ + (,s2 + ^ + 2 - 120 1.9^1' 



N 

-(6e-l)(.9.5s*+.9:.9s) } , (2.31) 



5 



and similarly for the T„„ component 



Ciu){T,,) = - ~ 40 + (1 - X - 8C + 4^x) p?] 

— 2 ~2 

+ (6C - + (1 - 20lp, + ^ + 



N 



oo 2 



2 ./o 27r2 



+ (6?-l)(5s.9*+g:5s) 



(2.32) 



It can be easily shown that the energy-momentum tensor is diagonal and that the spatial components are all equal, 
except for a geometrical factor (see Appendix . We then define [|| 



[ Tpjy ) = C(u)diag {e, p, p sinh^ rj, p sinh^ rj sin^ 6 ) . 



(2.33) 



III. RENORMALIZATION OF T^^ 

In order to analyze the divergences of we make use of the adiabatic expansion of the modes gs, since we know 

^ that (T^jy) and adiabatic t^^-ve the same ultraviolet behaviour. Up to second adiabatic order we have 



5s 



exp 



i / du'ns{u') 



(3.1a) 



1 

2n, 



1 

2wc 



8^4 



(3.1b) 



2n, 




16^3 



+ 



(3.1c) 



A - A* 



217? 



1 cO^ 
"2u;2 



(3.1d) 



where lo1{u) — + x{u). 

We regularize our integrals by introducing a non-covariant cutoff A in physical momentum, which corresponds to 
a comoving momentum cutoff Sm — C^^'^{u)K. We shall work with a fixed physical cutoff A, so that the comoving 
cutoff Sm depends on the conformal time u. 

We have three physical parameters to renormalize: the mass of the pion m.„, the self-interaction coupling constant A, 
and the (dimensionless) coupling constant to gravity ^ necessary in order for the field theory to be renormalizable, even 
in Ricci fiat spacetimes, such as the one considered in this paper p^ ]. The quadratic divergences in the gap equation 
will be removed by mass renormalization, and the logarithmic divergences will be subtracted by renormalization of 
the coupling constants A and ^. In addition, we shall see that there is an extra quartic divergence in both the energy 
density and the pressure, coming from the mode integrals in (T^i/) that carry an extra factor of fc^, which can be 
removed by renormalizing the cosmological constant. 

First we examine the divergences in the gap equation ( ^.28 ). From the adiabatic expansion (j3.l|), we know that the 
mode integral appearing in this equation has a quadratic divergence. In order to remove this divergence, we perform 
mass renormalization by subtracting the regularized gap equation for the vacuum state: 



27r2 



(s2 + c2«) 



,1/2 



(3.2) 
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This yields for the gap equation 



X{u) _ C{u)m^ 



A 



A 



+ pUu)-C{u)f^ + N 



s ds 
2^ 



(2n, + 1)15,^1' - 



2(s2 + e2«)i/2 



(3.3) 



Note that the second integral is independent of time. Once we have removed the quadratic divergences in the gap 
equation, we are only left with logarithmic divergences, which are subtracted by the following coupling constant 
renormalizations 



1 
A 



1 



N 
1 



1 



27r2 (s2 



(3.4a) 



Note that ^ = 1/6 = is a fixed point of the renormalization flow equations |12[ 
The renormalized form of the gap equation is then given by the following 



^ = -"«-^6^+''^+^y„ ^(2n. + l)|g.(.)| 

NC{u) k^dk 1 N[x{u) - C{u)m^] Pdk 



27r2 ^2)1/2 4 7o 27r2 (A;2+rn2)^/^ 



(3.5) 



We now proceed to write the regularized expressions for the energy and isotropic pressure and analyze the diver- 
gences appearing in the mode integrals. We shall take £,r = 1/& as suggested by [Q. The physics behind this choice 
is clear. If one considers an arbitrary composite scalar field (such as the meson field considered here) and an effective 
field theory at a scale A, and one carries out a renormalization group analysis in the leading large- approximation 
(or a fully improved one- loop renormalization group approximation), ^ = l/Q — ^^b. is found to be an attractive 
renormalization group fixed point in the infrared limit This means that even if there are corrections to ^ = 1/6 
at large scales (of order A), the observed low energy value of the coupling ^ tends to a physical value oi — 1/6, as 
one evolves into the infrared. The choice ^ = 1/6 = also implies that the divergences of both the energy and the 
pressure can be obtained at adiabatic order zero |p^ . 

Now we examine the divergences present in the mode integrals of the energy density and pressure. Recall that 

E{u) ^ C'{u)eiu) ^1[p1+ pUx " 1)] - Ip^ ^ ^ " ^ 



and 



2 ./n 27r2 



Pin) ^ C\u)p{u) EE 1 [p2 _ p2(^ ^ ^ 2.^^^ + ^ + ^ 

A^ s^ds 



{2n, + l){\gsf + {s'^x)\9sn- (3.7) 



6 Jo 27r2 

For the energy the contribution coming from the mode integrals is 

N s^ds 

~2 . L ~2 



^(2n,-t-l){|<7,|V(s2-fx)|5sp} , (3.8) 



and using the adiabatic expansion given in Eq. (3.1), the divergent part of this integral is given by 



A^ 



s'^ds 




s^ds 


2 '^'^ 




27r2 


Nslx 


Nx\ 
- 64.2 1- 




167r2 


K X 



Nst Nslx Nx^ fisl \ ^ (x' s 
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For the pressure the contribution coming from the mode integrals is 



TV 
~6 



271- 



-{2ns + l){\gs\' + {s'~x) \9s\'} 



(3.10) 



Again using the adiabatic mode expansion, the divergent part of this integral is given by 



C\u)po = ^ 



s ds 



1 



2^2 ^,2 



6 ^0 

487r2 487r2 ' Qin^ 



2 ~ Ar -2 



In 



X 



7Nx^ 
3847r2 



(3.11) 



Notice that eo ^ —po, as is required by general covariance (see Appendix We must then enforce covariance 
by hand. We mention here that the introduction of a momentum cutoff does not spoil covariance at the level of the 
energy density (see Appendix |^) . On the other hand, we must carefully handle the subtractions that will lead to the 
finite isotropic pressure [T3[ . The quartic subtraction in the energy density and pressure is a renormalization of the 
cosmological constant; in other words, a subtraction of the cosmological vacuum energy TVA'^/IGtt^. 

If we now define 



Nst 



167r2C2(u)' 



(3.12a) 



PR. = P - Pq - f a + 



Nst 



167r2C2(u) 



(3.12b) 



we can easily show that the energy density and pressure are finite, by making use of Eqs. (^j) and (| 

In order to numerically evaluate the physical energy density and pressure, we must also make sure that the vacuum 
energy that we are measuring with respect to is zero (i.e. when we are at the minimum of the potential). This vacuum 
energy is calculated in the "out" regime of the collision, when (a) = /rr, (7?^) = 0, x = 'J^Ji a-nd the mode functions 
are the zeroth order adiabatic ones with 0^2 = ^2^ C'{u)m^. The renormalized vacuum energy is given by 



2 f2 



m 
4A 



167r2 



N 

2C2,o 



s ds 



7^ / 7-^\/s2 +m'^C . 



2tt 



(3.13) 



The vacuum value of the pressure is Pvac — ^vac- These values must be subtracted from Eqs. ( |3.12[ ) so that when we 
have reached the "out" regime, the energy density and pressure go to zero. 

We calculate now the trace of the energy-momentum tensor, and show that these definitions of en and pn give a 
finite trace. We have from its definition 



If we make use of Eqs. (3.12), we can write 



NxHu) 
327r2 



where we have defined 



m'Ciu) Nm^C{u) 

+ LCm - 



Ns 



Xr 



167r2 



m ~2 



Nsl 



(3.14) 



(3.15) 



(3.16) 



The first two terms in Eq. (3.15) are the renormalized classical trace of the energy-momentum tensor, and the 
second term is the one-loop quantum trace anomaly. 
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IV. NUMERICAL RESULTS 



A. Energy-Momentum Tensor Conservation 

The renormalized energy-momentum tensor obeys the conservation law 

V^(Tnfl-0. (4.1) 
Using the ChristofFel symbols tabulated in Appendix we find that the v = u component of the conservation 



equation takes the form 



In terms of the variables 



eR{u)+i[pR{u) + eR{u)]=0. (4.2) 



Ns 



,4 



^B =£- -T^, (4.3a) 

Vn=V-C^{u){po + e^) + ^, (4.3b) 

we can rewrite the conservation equation as 

£r{u) + ?>Vr{u) -£r{u) = Q. (4.4) 



Recall that 



1 ~ \^ 

N f"^ S^ds ,^ . ,2 ,2 |2l 



Therefore 



2 ^^^(2n, + l){|g-r + (.^+x)|5.r}-Y^. (4.5) 



ffl = PiPl + (x - + -^XPt ^ -^^^ ~ 2A ~ ~^ ^'^'^ 

N s^ds 



2 ./o 27r 



{2ns + 1) [gsgl + gsijl + (s' + x)(g..9: + gsgl) + x\gs?] 



+ Y#(2"^™ + 1) + (^m - ^ , (4.6) 

where we have taken into account the contribution coming from the upper limit of the mode integral, since Sm depends 

on u. Making use of the identities Vi — Vi, — 3ji, the equations of motion, and using the adiabatic approximation 
for (we can always choose Sm to be a high enough comoving momentum so that the adiabatic approximation is 
valid in this limit), we then obtain 

£r = -X«' - 3j,p. + ^sl{2ns^ + l)ujs^ ~ ^ . (4.7) 



The occupation number n^^ goes to zero for large Sm, so we can neglect it. Expanding out yields 

. 24, 8^4 

So that we have 



= yMT^ = l + 7T^-^ + -- - - (4-8) 



87r2 327r 

or equivalently 



£r = -XV + ^JzPz - TJTT^ , (4.9) 



ffl = ~XVI - i]^P^ - 1^ . (4.10) 

If we compare this expression with Eq. ( [3.15 ), we can easily see that the conservation equation is satisfied. 
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B. Effective Equation of State 



Now that we have finite equations for the energy density and pressure, we can evolve them in proper time and 
numericaUy investigate their behavior. Since we have a nonequilibrium situation, it does not really make sense to 
calculate an actual equation of state, i.e. p = p{e). However, since e and p are both functions of u, in regions where 
both are monotonically increasing or decreasing, one can determine an "effective" equation of state p = p(e). In this 
calculation, we do not have two-body scattering, so that there are many oscillations in the pressure. We expect at the 
next order in the 1 /N expansion, when these effects are included that one will be able to extract an effective equation 
of state. 



C. Numerical Simulations 



We choose the initial state at a temperature above the phase transition in thermal equilibrium, with all particle 
masses positive. The equations are solved self-consistently at the starting time to obtain the values of the mean fields. 
We fixed the value of x the initial time as the solution of the gap equation in the initial thermal state. We also 
required that the initial expectation values of the a and tt fields satisfy 7f^(uo) -I-ct^(uo) = a"^, where ax is the thermal 
equilibrium value of $ at the initial temperature T. The critical temperature is Tc ~ 160 MeV, so we choose T = 
200 MeV, which gives ax = 0.3 fm The mode functions are chosen as their adiabatic values [see Eq. (2.24)]. We 
choose the coupling constant A_r — 7.3 [D, and take the conformal value of the coupling to gravity £, — ^n. — ^ / Q [0- 

In Fig. ^ we show the time evolution of the x field, which is the effective mass squared for the ip field. When 
this field becomes negative, long-wavelength unstable modes begin to grow exponentially. Therefore this field serves 
as a measure of instability in the system. For initial conditions with a negative derivative, we find that the x field 
is negative for about 3 fm/c of proper time. Figure |^ shows the numerical calculation of the energy density and 
pressure, in units of fm^"*, for initial conditions where no instabilities arose. Note that the system quickly reaches the 
"out" regime and relaxes to its vacuum values. Figure || shows the same plot for initial conditions with an instability. 
For the unstable initial conditions more energy density was produced initially, but both the energy density and the 
pressure drop off very quickly, as is expected (since the system is cooling very rapidly). It is interesting to note that 
the pressure becomes negative when x is negative. If we define the speed of sound Cq(t) by p = CQ(T)e, we notice that 
4<1. 



1.0 




0.0 1 ; 



-0.5 I I 

5 10 15 20 25 

T (fm) 

FIG. 1. Proper time evolution of the x field for the following initial conditions. Solid line; cr(ro) — ar, 7ri(ro) — 0, and 
^{tq) = 0. Dashed line: o-(ro) = ctt, 7ri(ro) — 0, and (T(ro) = —1 
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FIG. 2. Proper time evolution of the energy density and pressure for the initial conditions (T(ro) — (Jt, T^i{To) — 0, and 
a(ro) = 0. 




-1.0 



FIG. 3. Proper time evolution of the energy density and pressure for the initial conditions (T(ro) — or, 7ri(To) — 0, and 
a(ro) = -1. 

Figure || shows the conservation of the energy-momentum tensor. Notice that it is not conserved for short proper 
times. Because of the presence of the Landau pole it is not possible to take the physical cutoff A very large (we use 800 
MeV). Due to this rather small cutoff, the occupation number ns„^ is not actually ze ro, w hich it was assumed to be 
for the derivation of the energy- momentum conservation [see the discussion after Eq. (4^)]. As soon as the cutoff Sm 
becomes large enough, then the occupation number goes to zero and the energy-momentum tensor is conserved. Figure 
^ shows the energy-momentum tensor conservation for initial conditions with an instability. We see a qualitatively 
similar behavior, regardless of initial conditions. 



11 




-12.0 



10 15 
x(fm) 



25 



FIG. 4. Energy-momentum tensor conservation, for the initial conditions <7(ro) = ctt, 7rj(ro) = 0, and &{to) = 0. 
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FIG. 5. Energy-momentum tensor conservation, for the initial conditions cr(To) = ctt, ni{To) = 0, and &{to) = — 1. 
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V. CONCLUSIONS 

Wc have shown how to regularize and renormalize the energy-momentum tensor for a spherically symmetric model 
with a time-dependent comoving momentum cutoff which breaks covariance. We computed the finite energy density 
and pressure, and numerically examined the equation of state and proved energy conservation for this system. We 
found that the energy density and pressure both decrease quickly, which is expected for a rapidly cooling system. An 
interesting feature of the expansion is that when the effective mass squared is negative, the pressure is also negative. 
A negative pressure implies cavitation, which could mean the formation of domains, i.e. DCCs. 
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APPENDIX A: CHRISTOFFEL SYMBOLS 

The only non- vanishing Christoffel symbols for the metric 

9nv = C{u) diag (1,-1, — sinh^ 77, — sinh^ 77 sin^ 6) , (Al) 

C{u) = a2(r) =T^= e2"/m2 , (A2) 



with 



are 



= 1, r^, = sinh%, 

= sinh^ rj sin^ 6 , Tgg = — sinh 77 cosh 77 , 

= — sinh77C0sh77sin^ , F^^ = — sin^cos^, 



■pu _ pr/ P'P 1 

uu rju 9u i^u ' 

r^, = r^,, = coth77, 



= cot' 



The Christoffel symbols are used to derive the conservation law for the energy-momentum tensor. 



APPENDIX B: ENERGY-MOMENTUM TENSOR FOR A PERFECT FLUID 



For a perfect fluid we can write the energy-momentum tensor as follows 

T^ll, = -p g^v + (e + M^Wi. , (Bl) 

where are the components of the velocity vector of the fluid, such that it is normalized, i.e. g^^u^Ui, = 1, which 
corresponds to a timelike vector field. In the reference frame where the fluid is at rest we can write 

W/i = {Uu, Ur,, U0,U^) = {uu, 0) , (B2) 

and the normalization condition simply implies that g„„ = u,„7i„ = C{u). It is then straightforward to obtain the 
components of the energy-momentum tensor in this coordinate system 

Tuu = 9uu e = C{u) e, 

Ttiv = -9vv P = 

Tee = -gee p = C{u) psinh^ 77, 

T^p^p = -g^ip p = C{u) p sinh^ 77 sin^ 9 . (B3) 

Then we can write 

( T^i, ) = C (w)diag {e, p, p sinh^ 77, p sinh^ 77 sin^ 6 ) . (B4) 
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APPENDIX C: REGULARIZATION BY COVARIANT GEODESIC POINT-SPLITTING 



The regularization scheme (introducing a physical momentum cutoff A) used to render the mode integrals finite in 
the energy density and isotropic pressure is not a covariant method. Since we must obtain covariant results (we want 
to obtain (I)„y)ronormaiizcd), we havc to proceed with care, in order to perform the physical subtractions that will yield 
the physical finite energy density and pressure. We shall make use of the results obtained by Christensen Q|, in order 
to analyze the non-covariant structure on the divergences in eg and po, since we know that covariant point-splitting 
in a spatial direction is equivalent to regularization by introducing a momentum cutoff. 

Christensen B obtained: 



C^riLf) quadratic — liuij;/ — - o / \ \ I / \ \ 

An-' (ct^cta) 2 \ cta) 

{Ti_ii,) logarithmic = ( ^ 2 ' ^ 7 ' ^^^^ 



If we choose cr^ to be a spacelike vector such that 

1 



3 



gij{(T (7a) and cr„ = with i,j = rj,9,ip, (C2) 



we can write 



i'-^uu) quartic hm^^/ — 'X 9 / \ \2 ^(^)' 
^TT (ct-^CTa) 

(T ) -lim ^ ^ 

'jpi I quartic ^J^^-i-ia: 



"2^2(^A^^)2 3 ' 
J ^ ^ Cf. 

47r2 (a^CTA) 2 ^ 
1 1 X C(«) 



fix 

(T^uu) quadratic — liui^;/ — >x~. T7 \ 7 77 

477^ (ct^CTa) 2 

i'^rjr]} quadi^atic — hm^; 



4772 (ct^cta) 2 3 ' 

{Tuu) logarithmic = \ ^ ^ ^ ^ 

{Tnr,) logarithmic = ^ C {u) + ^ I ^ I ) ' (^2) 



^From Appendix ^ we know that T„„ = C{u)e and that T",,,, = C{u)p and from the previous equations (C3) we can 
make the following identifications: 

^quartic — ^l^a 



-V 1 1 1 

1 1 X 



'4772 (ct^CTa) 2' 

1 1x1 



^quadratic — lim^:' — 
Pquadratic — hui-j;' — 



^logarithmic — — Ihll^'^x ^p-j — ( 7 + | — 



4772 (ct^o-a) 2 3' 
4772 8 V ' 2 '2 



1 / 1 ^cr . 

Plogarithmic = YoTlx'^x^^ ~8~ V ^ 2 ' ~2~ ' ' ' 



Thus it follows 
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^quadratic 
^quartic 
^logarithmic 



3 Pquadratic: 
Pquartici 
Plogarithmic 



(C5) 



We can now compare these general results for spatial point-splitting with Eqs. (3^) and ( 3.11 ) by making use of 
the natural identification a^(J\ oc A"^. It is e asy to see that the quartic, quadratic and logarithmic terms of these 
expressions fulfill the requirements of Eq. (C5), as we wanted to show. 
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